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Abstract. We provide explicit formulas for the quantum integrals of a semi-infinite q-boson 
system with boundary interactions. These operators and their commutativity are deduced 
from the Fieri formulas for a q —>■ 0 Hall-Littlewood type degeneration of the Macdonald- 
Koornwinder polynomials. 
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1 Introduction 

The g-boson system [1, 24] constitutes an integrable g-deformed lattice regularization of the 
quantum nonlinear Schrodinger equation [11, 15, 17] built of g-oscillators [12, 20]. Its n-particle 
Bethe Ansatz eigenfunctions amount to the celebrated Hall-Littlewood functions [7, 16, 26]. The 
model in question can moreover be viewed as a degeneration of the recently found stochastic 
g-Hahn particle system [2, 21, 25]. 

By deforming the g-oscillator algebra at the boundary, a semi-infinite g-boson system was 
constructed [6, 8] with eigenfunctions given by the hyperoctahedral Hall-Littlewood func¬ 
tions [18, 27] that arise as a ((7 —)• 0) limit of the Macdonald-Koorwinder polynomials [14, 18]. 
In the present note, we use a corresponding q ^ Q degeneration of the Fieri formulas for the 
Macdonald-Koornwinder polynomials [3] to arrive at explicit formulas for the commuting quan¬ 
tum integrals of the latter semi-infinite q-boson system with boundary interactions. 

For the g-boson systems on the finite periodic lattice and on the (bi-)infinite lattice, analogous 
descriptions of the commuting quantum integrals stemming from the Fieri formulas for the 
Hall-Littlewood functions can be found in [5, 16, 26] and in [7], respectively. Freviously, Fieri 
formulas for Macdonald’s (g-deformed Hall-Littlewood) polynomials were interpreted in a similar 
vein as eigenvalue equations for the quantum integrals of lattice Ruijsenaars-Schneider type 
models [4, 10, 22, 23]. 

2 Commuting quantum integrals 

Let ei,...,en be the standard unit basis of and let A denote the cone of integer partitions 
A = (Ai,..., An) with parts Ai > • • • > An > 0. 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 
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Eor I G {1,... ,n}, we define the following difference operator Hi acting on the space C'(A) 
of complex lattice functions /: A —)• C: 


W)(\} := 5: Uj^nJt,i-\J+\-\J-\W^J+,J-Wfi^ + 6J+ “ ej-)) (1) 

J+nJ_=0, |J+|+|J_|<i 
A+e —e j_ SA 


where | J| refers to the cardinality of J C {1,..., n}, H := {1,..., n} \ J, and ej := 
coefficients of the difference operator in question are built of the factors 


E ej. 

j&J 


The 


1 _ J-l + fc— 

MG,._(A):=ni n (i-w”-A)) n 

\ 0 <r<s <2 . - 


Xj=Xf^ , Ej 


with ej = €j{J+, J-), and 

UK,mW := (-ir 




E 


1+,I-CK 


(l — (l — tot2t'^ 


i 6 /+ 


1 J^-L- 

l+nI-=0,\I+\ + \I-\=m \\j=' 

n 

jGl- 

Aj=l Xj — 


1 - 

1 - 


n 


Aj=A;c + l 


1 - 

1 - 


\ 


xfic’ n * n * ' 

j,kGK, j<k j,kGK,j<k . 

with ej = Cj(/+,/_). Here we have employed the notation 


£j('A|-, «/-) 


1 if j G J+, 

< -1 if j G J_, 
0 otherwise 

\ 


for J+, J- C {1,..., n} with J+ n J_ =0, and we have also assumed the standard convention 
that empty products are equal to 1. 

When ^ = 1, the action ol Hi (1) is relatively straightforward. Indeed, in this special case 
our operator amounts to a second-order difference operator of the form 

(ifi/)(A)=u(A)/(A)+ ^ ^+(A)/(A + e,)+ ^ wj(X)f {X - Cj), ( 2 ) 

A- 1 - 6 jGA a—gA 


with 


w+ (A) 


H (i-w^) 


L 0 <‘r<s <2 


n 

j<k<n 
Afc— Xj 


1 - 

1 - ’ 


(A) = 


n 

Xk=Xj 


1 - p-^ 


and 


u(A) 


l<j<n 
X-\-€.j G A 


j<k<n 

Xk—Xj 


1 _ ti-Hfc-i 

1 - 













Quantum Integrals for a Semi-infinite (;-Boson System 


3 


where 5^ •= 1 if = 0 and 6m ■= 0 if m ^ 0. Upon adding a harmless constant term of the 
form Y1 to the potential u{\), this reproduces the action of the semi- 

l<j<n 

infinite g-boson Hamiltonian from [ 6 , Section 4.1] restricted to the n-particle subspace (cf. also 
the proof of [ 6 , Proposition 3]). The g-boson Hamiltonian considered in [ 8 ] amounts in turn to 
the parameter degeneration t 2 —^ 0 (cf. [ 6 , Section 4.2]). In this interpretation the parts of A 
represent the positions of n interacting quantum particles - dubbed g-bosons - that hop on 
a one-dimensional half-lattice formed by the nonnegative integers. The parameter t corresponds 
to the deformation parameter q of the underlying g-oscillator algebra and the parameters t 2 

play the role of coupling constants for the boundary interaction on the lattice end-point at the 
origin (cf. Remark 2 at the end) [ 6 , 8 ]. Here we will generally think of these parameters as 
rational indeterminates, unless explicitly understood otherwise. 

The main result of this note is the following theorem, which states that the difference opera¬ 
tors Hi ,..., Hn constitute a system of commuting quantum integrals for the n-particle g-boson 
Hamiltonian H = Hi (2). 

Theorem 1. The difference operators Hi,, Hn (1) mutually commute. 


3 Proof 

Our proof of Theorem 1 hinges on the g —)• 0 degeneration of an explicit Fieri formula for the 
Macdonald-Koorwinder polynomials from [3, Section 6 ]. 

3.1 Macdonald—Koornwinder polynomials at qr = 0 

In the limit g —>■ 0, the Macdonald-Koornwinder polynomial [14, 19] gives rise to a two-parameter 
extension of Macdonald’s HC-type Hall-Littlewood function [18, § 10] of the form [27] 

PxiO = cxYl A G A, ( 3 ) 

w&W 

with 

3 

r iGi ■= rr (1 - (1 - yr 

i</i<n (1 - e'fe-^'^)) (l - l-e2*U 

Aj>0 

and 

^ TT __ 

l<r<3 

n 

Here (A,Q = Xjf^j and W = Sn ^ {— 1 , 1 }"' denotes the hyperoctahedral group of signed 
i=i 

permutations acting linearly on the components of ^ = (^i,... ,^n) G C". The trigonometric 
polynomial (3) is normalized so as to achieve a unit principal specialization value [ 6 ]: 

Px{ip) = f at p= ^ (log(to) + (n-j)log(t))ej. 

l<j<n 
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Eor parameter values in the domain 


G (- 1 , 1 ) \ { 0 }, r = 0 , 


the g —)• 0 degeneration of the Macdonald-Koornwinder orthogonality ensures that [14, 27] 



Px{op,iom)\"d^ = o 


if A / /i, 


where the orthogonality density is given by the squared modulus of 


A (0 


n 

l<j<k<n 


(l — (l — 


n 


1 — 


l<i<n Yl (l-trpij) 

r=0 


(4) 


3.2 Fieri formulas 

Eor a given choice of generators Ei{^),... ,En{C) for the algebra of trigonometric polynomials 
with hyperoctahedral symmetry, the associated expansions of products of the form Ei{^)Px{^) 
in the basis P^{^), /U G A give rise to a system of recurrence relations commonly referred to as 
Fieri formulas. 

Theorem 2 . For I G {l,...,n}, the hyperoctahedral Hall-Littlewood function Px{(,) (3) with 
fa = 0 satisfies the following Pieri-type recurrence relation: 

Eimxd) = E UjlnJt,l-\J+\-\J-\Wyj+,J-WP\+ej^-ej_ iC), 

J+,J-C{l,...,n} 

J+nJ_= 0 , |J+|+|J_|<z 
A+ej^-ej_eA 


with 


l<ji<-"<jl<n l<k<l 




ieJ+ jeJ- 

A,=0 A,=l 


X 


n 

l<_ 7 <fc<n. 
Xj =Xj^ ,6j 


1 - 


l<_ 7 <n 


where ej = ej(J+, J-), and with UK,m{^) being defined as in Section 2. 

Proof. The stated formula boils down to a degeneration of the Fieri formula for the Macdonald- 
Koornwinder polynomials in [3, Theorem 6.1]. This degenerate formula is obtained through 
a straightforward but somewhat tedious computation that involves setting t^ = q and performing 
the limit g —)■ 0. The present formulation moreover employs a compact expression for the 
mnltiplying polynomial E'/(^) stemming from [13, equation (5.1)] (cf. also [9, Section 2]). ■ 

When t 2 = 0, Theorem 2 reduces to a Fieri formula for Macdonald’s BC-type Hall-Littlewood 
functions. 
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3.3 Commutativity 

For A G A such that A -|- ejj_ — ej_ G A, we have the following functional relation between 
the coefficients of the difference operators in Theorem 1 and those of the Fieri formulas in 
Theorem 2: 

Vj^,j_{X)h{X + ej^-ej_) = Wj^,j_{X)h{X), 


with 


h{X):= 

l<j<n 

A>o 


The upshot is that upon conjugating Hi (1) with the (invertible) multiplication operator in C'(A) 
of the form / —)■ hf, the coefficients IFj^_j_(A) get replaced by the Fieri coefficients Vj^^j_{X). 
The Fieri formula in Theorem 2 tells us that the resulting conjugated difference operators 
commute on the joint eigenbasis of hyperoctahedral Hall-Littlewood functions (3) (viewed as 
lattice functions of A G A depending on a polynomial spectral parameter ^ G M"). Indeed, 
from this perspective the Fieri formula corresponds to an eigenvalue equation with the bounded 
function Ei(^) playing the role of the eigenvalue. The orthogonality relations (4) moreover 
guarantee the completeness of these (generalized) eigenfunctions (cf. Remark 1 below), i.e., the 
difference operators in question commute in fact as bounded operators in the Hilbert space 
^^(A, u) C C'(A) determined by a discrete measure with weights 

^x=([ |PA(OA(OPdC^ , AgA. 

W[ 0 , 27 r]" / 


This means that the operators commute in particular on the (stable) subspace of C'(A) consisting 
of the lattice functions with finite support. But then the commutativity must actually hold 
on the whole space C'(A), as given any / G C'(A) and any A G A, the evaluation at A of 
the commutator of two of such difference operators acting on / depends manifestly only on 
evaluations of / at a finite number of lattice points in A. Finally, the commutativity is extended 
beyond the parameter values in the orthogonality domain by analyticity. 

Remark 1. Notice that away from the boundary (i.e., for Ai > • • • > An > 0) the wave 
function P\{^) (3) decomposes as a linear combination of plane waves (of Bethe Ansatz form). 
In particular, the wave function in question does not belong to the Hilbert space and 

constitutes in fact a generalized eigenfunction of the discrete difference operators arising from 
the Fieri formulas. (The spectra of these bounded difference operators are absolutely continuous 
rather than discrete.) The orthogonality relations (4) do nevertheless imply that any / in u) 
can be represented through a wave packet via the associated (generalized) Fourier transform: 

/(A)=/ AgA, 


where 


AeA 


(5) 


which confirms the completeness of our generalized eigenfunctions. (Here the convergence of 
the sum on the r.h.s. of (5) is in the strong L^([0, 27r]"', |A(^)pd^) Hilbert space sense.) 
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4 Extension to four-parameter boundary interactions 

In principle there is no genuine obstruction preventing us from adapting the commuting quan¬ 
tum integrals of Theorem 1 to the case of the more general semi-infinite g-boson system in [6] 
with four-parameter boundary interactions. Eor this purpose, one needs to establish a ge¬ 
neralization of the Fieri formula in Theorem 2 covering the hyperoctahedral Hall-Littlewood 
polynomials P\{^) (3) with arbitrary. It turns out, however, that in this more general set¬ 
ting the coefficients of the Fieri formulas (and thus also those of the corresponding quantum 
integrals) become quite baroque. We wrap up by indicating briefly how the above formulas are 
to be modified when dealing with such general boundary interactions involving four parameters 
to,. . . , ts- 


4.1 Fieri coefficients 


Even though the global structure of the Fieri formula for the hyperoctahedral Hall-Littlewood 
functions Px{^) (3) remains of the form described by Theorem 2 when dropping the condition 
that ts be zero, the fine structure of the coefficients is now more intricate: 

(l - rt^-j+moW+m^{\)-m+{X)'^ 0 (1 “ 

yp,j-w=ll - 

i 6 J+ 


(1 - rt2(«-i))(l - Tp{n-j)+l) 

(I - n (1 - 


(1 


l<r<s<3 


i 6 J+ 


ieJ- 


(I - rt2{’^-i))(l - 


1 fl+k—j 

n n 

l<j<k<n 
Xj =A/(., Ej 


and 


= (- 1 )’ 


E 


I+,I-CK 
I+nI-= 0 ,\I+\+\I-\=m \Aj=i 


n (1 - tout’'-’) 


jei+ 


i 6 /+ 

Aj=l 

I — ^fin+l-j-k 


n l<r<s<3 T—r 

li l_^i2n-j-k 

j£l- j,kGK,j<k 

Xj=l ej-|-efcG{—2,1,2} 

Aj=l) Xk—^l + ej^ 


-|-r 1 - 

11 1 _ ^f2n-j-k 11 I _ fk-j 11 I _ fk-j 

jGJ+U/_,/cGX\7_ jGJ_,A:G/+ 

j<.k, Cfc—ejG{0,l} Xj=Xk-, ej>efc Aj=Afc-)-l 

Aj=(5i_^e^., Afc=0 


n ‘P n *n ‘' 

jGK j^k^K, j<k j^kGK,j<k . 


( 6 ) 


where mi{\) = |{1 < j < n | Aj = Z}|, m^(A) = |{j G J+ | Xj = 01) r = These 

Fieri coefficients are obtained from [3, Theorem 6.1] in the limit g —0. 
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4.2 qr-Boson quantum integrals 

By virtue of the argumentation in Section 3.3, the Fieri formulas at issue give rise to commuting 
difference operators Hi ,..., Hn of the form stated in equation (1) with 


j_ (A) = n (1 - H 




l<j<k<n 
Xj =A/», €j 


1 - 

1 - t^-i 


n 

i 6 J+ 


0<r<s<3 

(1 - rt2{^-i)-l)(l - rt2(’T’-i))2(l _ Tf2{n-j)+l^^ 


and with UK,m{^) taken from equation (6). The relation between the coefficients of the Fieri 
formula and those of the difference operators is again governed by a functional identity of the 
type in Section 3.3 with 

h{X)= Yl {I - . 

l<r<s<3 

When I = 1 the corresponding difference operator Hi now reduces to Hi (2) with 

1 - 


w 


+ (A) = (l - 


j<k<n 

X}g=Xj 


1 - 


0<'r<s<3 


(1 - - Tt2(»^-i))2(l _ 7-i2(n-j)+l) 


W, 


^ 1 _ A+3-k 

■W= n -rriFV’ 

l<k<j 

X]^=Xj 


(7a) 

(7b) 


and 


n(A) = - 


l<j<n 
X-\-G'j gA. 


/ n (i-toM’^-Q 

l<r<3 


5a, 


(1 - rt2(’^-i))(l - Tt2("-i)+i) 


n 

j<k<n 

Xi^=Xj 


1 - 

1 - t'^-3 




Y tot' 


n-j 


/ (1 - Tt^-^- 1 ) n (1 - trtst^-^) \ 

l<r<s<3 


(7c) 


A—e -j' G A 


V 


(1 - rt 2 (’^-i)-i)(l - Tt 2 (»^-j)) 


1 i 1 _ y-fc 

l</c<_7 

Xk=Xj 


Upon adding the constant term Y1 {tot^ •^)), this reproduces the n-particle g-boson 

l<_7<n 

Hamiltonian of [6, Section 3]. 

Remark 2. In [6] a Fock space description of the particle Hamiltonian Hi (2), (7a)-(7c) was 
provided as a system of q-bosons on the nonnegative integer lattice perturbed at the lattice-end. 
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Specifically, the boundary interactions arise in this picture from a deformation of the g-boson 
field algebra at the origin and its nearest neighboring lattice point parametrized by to, ti, t 2 , to. 
In general, i.e., when r = / 0, the deformed g-boson field algebra is no longer ultralocal 

at the boundary as the commutativity between the creation and annihilation operators at the 
origin and its nearest neighboring site is lost. When to = 0, the deformation of the g-boson 
field algebra is restricted only to the lattice end-point at the origin and the ultralocality is 
restored [6, Section 4.1], When both t 2 = to = 0, the boundary interaction degenerates and 
decomposes into an interaction arising from a one-parameter deformation of the g-boson field 
algebra and a one-parameter additive potential term of the Hamiltonian, supported at the lattice 
end-point [8]. 
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